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O • Abstract 

We show that the ratio of different structure functions have a maximum which depends 
on energy. We argue, using the Golec-Biernat and Wusthoff model as well as the eikonal 
approach, that these maxima are functions of the saturation scale. We analyze leading and 
higher twist contributions for different observables to check whether a kinematic region 
^ \ exists where high parton density effects can be detected experimentally. 
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1 Introduction 



In this letter we continue our study of phenomena associated with gluon saturation in deep 
inelastic scattering (DIS) following our original study in Ref. [|IJ. Our goal is to find the 
kinematic region where the evolution of parton densities can no longer be correctly described 
by linear evolution equations^, |3|, ^, H]. To solve this problem analytically one has to find a 
solution of the non-linear evolution equation which was derived in Refs. |], [7|][] and in Ref. || in 
different approaches. In spite of the fact that the asymptotic solutions are known ||, a detailed 
description of the boundary of the saturation region is still lacking. We would like to determine 
how gluon saturation can be detected in current and future experimental data. 

The idea suggested in Ref. []J was that there is a breakdown of the Operator Product 
Expansion as one gets sufficiently close to the saturation region. This enabled us to find the 
scale Q 2 where the twist expansion breaks down. Since, at sufficiently large Q 2 , the saturation 
scale is the only scale of the problem, it should be expressed unambiguously in terms of Q 2 . 
We refer the reader to our previous publication for more details. 

In Ref. [[|] we used the eikonal model 0, [II], [12|, [L3|, |L4| to calculate the various structure 
functions relevant in DIS on nuclei, and then to expand them into the twist series. In this letter 
we will address the subject of DIS on a nucleon and apply the same approach developed for 
DIS on nuclei. 

The paper is organized as follows: We start with a brief review of our approach without 
going into details of the calculation (Sec. 2). In Sec. 3, we discuss numerical estimates regarding 
the nucleon target. Finally, in Sec. 4 we conclude by proposing possible ways to measure the 
different twists at HERA and LHC 



2 The Model 



As was pointed out in the introduction we use the eikonal (or Mueller-Glauber) model to 
calculate shadowing corrections. A complete description of the eikonal approach to shadowing 
corrections (SC) was given in Ref. |L4] . 



It was shown in Ref. [jlj, [15], |16[ that the eikonal 
approximation leads to a reasonable theoretical as well as experimental approximation in the 
region of not too small xb- There are weighty arguments in favor of using the eikonal model 
to calculate contributions of different twists [|IJ in various cross sections. This simple eikonal 
model is a reasonable starting point for understanding how the shadowing corrections occur. 
The use of the Glauber-Mueller approach was justified for a nuclear target 0, and it should 



hold for a nucleon as wellfl7[ 

Following Ref. Q and assuming that all dipoles in the virtual photon wave function have the 
size r\ = 4/Q 2 and neglecting the quark masses we obtain for the longitudinal and transverse 
photon cross sections 
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1 Strictly speaking this equation is derived for a nucleus target, however there are reasons to assume [[T^ that 
it holds for the nucleon as well. 
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where 3 Fi is the generalized hypergeometric function, G\\ is the Meijer function (see Ref . [p"8] ) , 
and we have defined^ 

ra ( 2 , _2na s x B G{x Bl Q 2 ) . . 

y S { X B,V ) = ^ ' ^ 

where G(xb,Q 2 ) denotes the gluon distribution in the nucleon, R the radius of the nucleon, 
and Q 2 is used as a shorthand for Q 2 {xb, Q 2 ) defined in Eq. ([|). Q 2 is a saturation scale in the 
framework of our model. 

Diffractive cross sections can be calculated using the simple relations following from the 
unitarity constraint: 

°(i*p)?,M/Q 2 ) = <i*p)tAQVQ 2 ) - 1 -v{i*p)tAzQ 2 JQ 2 ) (4) 

Formulae Eq. (JTj) , Eq. (g) and Eq. (f|) make it possible to estimate the contribution of any 
order twist to the particular cross section by closing the contour of integration in a Meijer 
function over the corresponding pole. 

We argued in the Ref. QTJ] , that studying ratios of the cross sections can yield very useful 
information about the saturation region. While the cross sections receive large corrections 
from the next-to- leading order of perturbative expansion, their ratios do not[[l4|, In this 
paper we will consider two ratios (Jl/ct and a® /cr®- Both of these ratios display a remarkable 
property as a function of photon virtuality Q 2 . At small Q 2 the ratios vanish since the real 
photon is transverse. At large Q 2 they vanish as well, as implied by pQCD. This leads us to 
suggest that both ratios have maximum Q 2 max - We would expect that this maximum is a certain 
function of the saturation scale, since at small xb and large Q 2 this is the only scale important 
in the process of dipole scattering. So, studying the behaviour of the Q 2 max with energy we can 
determine the Q 2 (xb) as well. Two problems occur. First, the confinement scale affects the 
position of the maximum of ratios. We can only hope that at large energies (small xb) the 
maximum will occur at sufficiently large Q 2 where this influence is small, though as we will see 
below, it cannot be completely excluded. 

Second, we do not know the precise form of the function Q 2 = f{Q 2 na ^) which has to be 
determined using the full twist expansion. We can, however, roughly estimate the form of 
this function using a simple approach. Namely, we illustrate the most important features (for 
finding the Q^ax) °f Eq. (ffh and Eq. (|j) by the following formulae which have the correct 
asymptotic behaviour for <c m 2 < Q 2 and Q 2 ^> Q 2 > m 2 , and which are valid for any 
model: 



Q 2 + Q 2 
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Q 2 + m) Q 2 + Q 2 
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Q 2 + Q 2 ' 
Q 2 Ql 



Q 2 + m 2 (Q 2 + Q 2 ) 2 



where the logarithmic contributions are neglected. Note, that we explicitly assumed here that 
at large energies the only relevant scale of the process is Q 2 . 



2 In our previous work Ql) we used different notation AC = Q 2 S . 
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The ratios of the cross sections then read 



R 



a T Q 2 + mj 



Q 2 Q 2 



Q 2 + m 2 Q 2 + Q 2 ' 

Taking the derivative of R and R D with respect to Q 2 and equating it to zero, we find that the 
ratio R(Q 2 ) has no maximum, while the maximum of the ratio R D (Q 2 ) is at Q 2 max = JmjQ 2 . 

There is a maximum of R{Q 2 ) which appears when we take the logarithmic contributions into 
account (R decreases logarithmically with Q 2 ) and turns out to be proportional to Q 2 . Our 
simple arguments show that 

1 . The maximum of R is much shallower than that of R D 

2. The maximum of R D is proportional to \J~Ql ■ 

3. Existence of the non-perturbative scale m 2 is crucial for appearance of maxima. 

The maxima that we have just discussed in general occur due to the well-known limiting 
behaviour of the longitudinal and transverse cross sections at small and large Q 2 . What we 
showed is that at small enough xb those maxima are situated at a scale, which can be expressed 
unambiguously in terms of the saturation scale. 

To check these conclusions, we plot the Q 2 max and Q 2 versus logl/xe in Fig. 0(a). Recall, 
that we obtained formulae Eq. ([!]) and Eq. @ by assuming that all dipoles in the virtual photon 
wave function are of the size 4/Q 2 . Thus, to get the correct (up to logarithmic contributions) 
numerical value of the saturation scale at fixed we evaluated the gluon distribution at 
sufficiently large value of Q 2 (in fact we took Q 2 = 100 GeV 2 ) where to a good approximation 
it is a constant. 

Note, that in making the simple estimations above we did not rely on a particular model. 
So, we also check our conclusions using the Golec-Biernat-Wusthoff model fl9| which is similar 



to ours, but for which we know exactly what the saturation scale is. Using their a we calculated 
the ratios R and R D and found the maxima. The result is shown in Fig. 0(b). Surprisingly, 
both results agree quite well with our rough estimations^] (estimations of Fig. 0(a) are valid at 
not too large logQ 2 )- Moreover, the maximum of R is proportional to Q 2 in contrast to the 
maximum of R D . This result can also be obtained from our simple estimation by including 
logarithmic corrections in Ot- 



3 Numerical estimations 

We now present the results of the numerical calculations. We performed the calculation using 
GRV'94 parameterization pT|, [22| for the gluon structure function xbG(xb, Q 2 ) for two reasons. 



First, we hope that the non- linear corrections to the DGLAP evolution equation are not included 
in it. Second, it enables us to perform calculations at values of Q 2 as small as 0.4 GeV 2 . 

We noted in Ref. |l[ that the serious drawback of GRV'94 is the existence of a domain in 
Q 2 where it gives the anomalous dimension of xbG(xb, Q 2 ) larger than 1. To overcome this we 
use the following amended gluon structure function 

x B G(x B ,Q 2 ) = x b G grv (xb,Q 2 )9(Q 2 -Q 2 (x b )) 

+ *bG grv (x b , Q 2 {x b )) 9(Q 2 (x b ) - Q 2 ) (5) 
Q 2 {xb) 



3 Recently it was shown in Ref. [g(J that the experimental data at small xb indeed exhibit scaling, the total 
cross section depends on the unique scaling variable r = Q 2 /Ql- 
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where 9 is a step function, and Q 2 {x B ) is a scale at which the anomalous dimension equals 
1 (it is a saturation scale within our model accuracy). Eq. (|5D can be used as a simplified 
parameterization of the behavior of the gluon structure function in the saturation region. Our 
results turn out not to be sensitive to the precise form of the behaviour of xG(x B ,Q 2 ) in the 
vicinity of the saturation scale. 

In order to compare our calculations with the experimental data one has to multiply the 
theoretical formulae for the cross sections (structure functions) by factor » 0.5. This factor 
stems from the next-to-leading order corrections which can be taken into account by modeling 



the anomalous dimension as described in Ref. 

To estimate contributions of different twists to various structure functions defined by 



F(x B ,Q 2 ) = -^a( 1 *p) 



we performed calculations using Eq. ([[]), Eq. (fj) and Eq. ([|). Note, that Eq. ([[]) and Eq. (|2]) 
are valid only in the massless limit m 2 =0. As was explained above, this is a reasonable 
assumption as long as we are not concerned with the behavior at small Q 2 . All calculations 
are made for the THERA kinematical region. The results of calculations for different twist 
contributions are shown in Fig. [I]. Similarly to the nuclear target case the following important 
effects are seen, when we compare contributions of the first two non-vanishing twists: 

1. In all structure functions except Fl there is a scale Q 2 for the particular structure 
function where these two twist corrections are equal. Q 2 is largest for F®; 

2. Q 2 grows as ib decreases; 

3. Twist-4 corrections in Fl cancels numerically with twist-4 in[] Ft in the region where 
their separate contributions are quite large. This is a reason to measure the polarized 
structure functions instead of the total one. 

In analogy with the scale Q 2 max the scale at which OPE breaks down is expected to be a 
certain function of the saturation scale, this function differs for the various structure functions. 
The question of what the analytic expression of this function is warrants further study. We, 
however, can assert that the scale at which the OPE breaks down is not smaller than Q 2 . 
The latter can easily be estimated by equating the analytic expressions for twist-2 and twist-4 
corrections obtained in Ref. The result is Q 2 = f ■ Q 2 up to logarithmic corrections, where 
the numerical factor / is of order unity (e.g. for ctt, / = 0.79). 

In Figs, |2| we plot the ratios F L /F T (x Bl Q 2 ) and FjP / F® (x B , Q 2 ) versus Q 2 . The calcu- 
lations were performed without neglecting the quark mass m/. As is expected on theoretical 
grounds the ratio has maximum which increases as ib decreases 0. In the kinematical region 
of THERA it is expected, that Q 2 ~ 1/xb[[J. The maxima in Fig. |2] occur at Q 2 max 3> m 2 
which implies that we are sufficiently far from the confinement region. The maximum of the 
ratio of the diffractive structure functions R D is more obvious, as we saw in sec. 2. 

We show in Figs. 0(a) the behavior of Q 2 nax {xB) as a function of x B . The increase is still 
smaller than predicted theoretically. Our conclusion is that the kinematical region of THERA 
is intermediate between the linear and saturation regimes. 



4 Conclusions 

The aim of our work is to understand to what extent the high parton density regime of QCD 
could be detected in DIS experiments at HERA, there the density of partons could be sufficiently 



This effect was also noted in Ref. [23 where another, though similar, modclpjj] for SC was used. 
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large to make SC significant. Using the eikonal model we observed two manifestations of the 
saturation scale: 

1. The position of the maxima Q 2 max in ratios F L /F T and F^fF® is intimately related to 
the saturation scale Q 2 {xb)- Q 2 ma x shows xb dependence typical for the saturation scale 
(see Fig. 0). On the contrary, the ratios of structure functions satisfying DGLAP equation 
have no maxima, as is readily seen from the leading twist formulae given in Ref . [p]] . 

2. The maximum of the ratio F^fF® is much more apparent than that of Fl/Ft, since 
the former occurs due to leading twist behaviour of the structure functions, whereas the 
later is due to logarithmic corrections to this behaviour. 

3. The results of our calculations show that there exists a scale at which the twist ex- 
pansion breaks down since all twist corrections become of the same order. We can use 
the DGLAP evolution equations only for Q 2 larger than this scale. The energy (xb) 
dependence of this scale suggests that the experimental data for deep inelastic scattering 
on a proton will help us separate leading and higher twist contributions. In the case of 
nucleon deep inelastic scattering such a separation appears to be a rather difficult task 
and has not yet been performed. 
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Fi gure 1: Different twist contributions to the various structure functions for DIS on the proton: leading twist 
(at high Q 2 ) - dashed line, next-to-leading - dotted one, exact structure function - solid curve. 
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